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Pion interferometry results from pp and dAu collisions at yfs = 200 GeV are presented using data 
from the STAR detector at RHIC. Radii from three-dimensional correlation functions are studied 
as a function of transverse mass. The HBT results from pp and dAu collisions are compared to pion 
HBT from AuAu reactions, also taken at STAR. 



Hanbury-Brown Twiss interferometry has long been 
used as a tool to characterize the space-time freeze-out 
geometry of particles produced in elementary particle 
and heavy ion collisions EL El HQ- At RHIC (Relativis- 
ts Heavy Ion Collider), there is a unique opportunity 
to study pion HBT in three separate systems, pp, dAu, 
and AuAu, using the same detector at the same center 
of mass energy. 

In this work, presented as a poster at Quark Matter 
2004, pion HBT results from pp and dAu collisions are 
compared to AuAu results all taken at yfs = 200 GeV 
using the STAR (Solenoidal Tracker at RHIC) detector. 
Of primary interest in this report is to compare how the 
three-dimensional regions of homogeneity vary with the 
transverse mass of the pion pairs in each system (Mr = 
y^m 2 + krp where kr is the pair's transverse momentum). 

The central observable in HBT is the two particle cor- 
relation function, C2. This correlation function is sen- 
sitive to several interesting physics effects: the quan- 
tum statistics obeyed by the pairs being studied (Boson, 
Fermion, etc.), the quantum field configuration (thermal, 
coherent, etc.), the source geometry, the source dynamics 
(space momentum correlations, flow, jets, etc.), and pair- 
wise interactions (Coulomb, strong, etc.). The broader 
goal of HBT studies in hadronic physics is to isolate and 
understand these individual effects in various systems un- 
der different conditions. 

Experimentally, the two-pion correlation function is 
obtained from the ratio C2(q) — A(q)/B(q). The quan- 
tity A(q) is the two-pion distribution of pair momentum 
differences, q = pi — P2- The mixed background, B(q), 
is the two-pion distribution generated from pion pairs 
coming from different events. 

The correlation function was measured in three dimen- 
sions using the Bertsch-Pratt decomposition of q: g out , 
q S ide, and gi on g- Here, q\ ong is parallel to the beam. The 
component q s i^ e is perpendicular to both q\ ong and to 
the transverse momentum vector of the pair, kx- The 
third component, # ut> is perpendicular to both <2i on g and 
g S ide- The measured correlation function was then fit to 
a Gaussian 

C2(q) = 1 + Ae _ ^out^out - ^siderfide - -ftfongrfong (l) 

where A represents the strength of the correlation at 



zero momentum difference. It is also commonly known 
as the "chaoticity parameter." The correlation func- 
tion was also studied in one-dimension versus Qinv = 
\/\(Ei — E2) 2 — (pi — P2) 2 1 • Here, the correlation func- 
tion was parameterized with the Gaussian 

C 2 (Q inv ) = l + Ae- ii2< 3-. ( 2 ) 

The primary source of systematic error in the pp and 
dAu systems comes from the assumption about the shape 
of the correlation function. Future work will include 
extensively characterizing the correlation function using 
other methods such as the Edgeworth expansion, Legen- 
dre expansion 5], and imaging 6]. Systematic error on 
the HBT parameters are estimated to be as large as 20% 
using a Gaussian shape assumption. In spite of this, the 
Gaussian shape still reasonably characterizes the global 
properties of the correlation function for purposes of com- 
parison between the systems. 

The primary subdetector used at STAR in this analysis 
is the main Time Projection Chamber (TPC), described 
elsewhere 7]. Particle identification (PID) was achieved 
by correlating the magnetic rigidity of a track with its 
specific ionization [dE/dx). For the pp and dAu analy- 
ses, a one sigma cut was made around the pion Bethe- 
Bloch band while the electrons, kaons, and protons were 
suppressed at the two-sigma level providing a pure sam- 
ple of charged pions. The contamination due to electrons 
and kaons is small and contributes to a slight systematic 
error of about 1%, an estimate based on changing the 
width of the sigma cuts. Good pion PID was achieved to 
a momentum of about 800 MeV and the TPC itself im- 
poses an lower momentum cut at 100 MeV. Pion tracks 
were selected at mid-rapidity, \y\ < 0.5, for this analysis. 

Before cuts, approximately 10 million minimum bias 
pp events and 8 million minimum bias dAu events were 
used. 

Two sources of pair wise systematic error, track split- 
ting and track merging, were addressed. Track splitting 
occurs when a single track is reconstructed as two in the 
TPC and track merging is when two tracks are mistaken 
as one. Both effects can affect the shape of the correlation 
function at very low q. In AuAu collisions, the correlation 
function is quite narrow because the pion source size is 
large (~ 6 fm). For that system, this makes eliminating 
splitting and merging particularly important. However, 
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in dAu and pp collisions, the source is comparatively 
small, giving a wide correlation function. In pp collisions 
at STAR, splitting and merging effects only significantly 
affected the first 30 MeV/c of a wide correlation function 
(~ 200 MeV/c). However, for future studies where a de- 
tailed understanding of low-q behavior will be explored, 
removal of these systematic errors will become increas- 
ingly important and so were implemented in this study 
for both the pp and dAu systems for completeness. 

A standard pairwise Coulomb correction was applied, 
similar to that done previously at STAR |g. This stan- 
dard Coulomb correction has been shown to over-correct 
the correlation function in heavy ion physics and is inad- 
equate in that context. However, the overall effect of this 
standard Coulomb correction is small for small sources. 
It typically affects only the first 40 MeV/c of the cor- 
relation function. The use of other correction schemes 
does not significantly affect the trends of the extracted 
Gaussian parameters for pp and dAu collisions. However, 
for future detailed shape studies, the Coulomb correction 
will become extremely important. Ideally, one would like 
to move away from implementing Coulomb corrections 
altogether and employ methods such as imaging jy| that 
incorporate all appropriate final state interactions into 
the extraction of the source function itself. 

The Bowler- Sinyukov- CERES scheme is currently the 
Coulomb correction method of choice for heavy ion HBT 
at STAR 8]. However, the appropriateness of this 
scheme in pp and dAu collisions is still being studied. 
For example, the Bowler scheme assumes that the source 
is perfectly chaotic and any deviation of the chaoticity 
parameter, A, from unity is the result of contamination, 
long lived resonances, or any other undesirable effects 
that artificially reduce the correlation strength. However, 
this assumption may not be appropriate for pp collisions 
where coherence and core-halo considerations could be 
important. 

Figure shows a one-dimensional tttt correlation 
function for a pp collision at STAR with dN/dr] = 4, 
representing a typical event. The darker data points, as 
indicated on the figure, are Coulomb corrected. Both a 
Gaussian and an exponential were fit the to the Coulomb 
corrected data. The corrected correlation function is 
clearly not Gaussian. As a guideline, the width of the 
correlation function in relative momentum space is in- 
versely related to source size. For an incoherent source 
of identical, spinless bosons, C^q — 0) should approach 
2. This in turn should be reflected in the parameter A, 
which should approach unity for a fit functions of the 
form of Eq. (0. Because of experimental effects (such 
as PID contamination) and physics effects (such as co- 
herence and resonances) the value of C2(q = 0) is often 
measured to be less than 2. 

It is interesting to compare C2(q) between different 
systems. Figure shows three tttt correlation functions 
versus Qi nv for pp (minimum bias), dAu (minimum bias), 
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FIG. 1: A typical one-dimensional tttt correlation function 
from a pp collision at s/s = 200 GeV/c 2 at STAR. Both Gaus- 
sian and exponential fits to the Coulomb corrected data are 
shown. The darker points are the Coulomb corrected data. 




FIG. 2: Typical one-dimensional tttt correlation functions 
for the pp (upper), dAu (middle), and AuAu (lower) systems 
at STAR. The pp and dAu correlation functions are from 
minimum bias data while the correlation function from AuAu 
is shown for central collisions. All are subject to the same kr 
cut. 



and AuAu (central) that have not been Coulomb cor- 
rected. As the system size decreases from AuAu to pp, 
the width of the correlation function increases. Note that 
Bose-Einstein correlations are not a small effect in dAu 
or pp collisions in comparison to AuAu. 

Figure © shows one-dimensional projections of 
a three-dimensional Bertsch-Pratt correlation function 
along the out, side, and long directions for pp collisions. 
The projections are 80 MeV/c wide in the "other" direc- 
tions for each plot. The hole, as seen in the out direction, 
is a kinematic effect. It is the result of a pairwise vector 
kr cut (0.15 < kr < 0.25 MeV/c) in combination with 
an implicit single particle pr cut (pr > 0.1 MeV) due to 
acceptance in the TPC. 

The Mt dependence of the three dimensional HBT 
parameters provides information about the space- 
momentum correlations present in each system. A trans- 
verse momentum dependence of the longitudinal radius 
along the beam axis is associated with a boost-invariant 
expansion. Such behavior is commonly seen in a variety 
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FIG. 3: Typical out, side, and long projections of the 
three-dimensional tttt correlation functions for pp collisions at 
STAR. The one-dimensional projections shown are 80 MeV/c 
in the other directions. The origin of the hole seen in the out 
direction is discussed in the text. 
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FIG. 4: The transverse mass dependence of the 3D Bertsch- 
Pratt radii in pp collisions at STAR. The line guides the eye 
with a 1/V Mt parameterization. 



of high energy physics experiments where HBT parame- 
ters are studied over a boost-invariant region of longitu- 
dinal phase space. Other space-momentum correlations, 
such as longitudinal flow in heavy ion collisions, become 
augmentations to this intrinsic dependence. A transverse 
momentum dependence of the transverse HBT radii is 
indicative of an explosive expansion of the pion source, 
usually associated with global collective properties such 
as transverse flow. In elementary particle collisions, a 
transverse momentum dependence of the transverse radii 
have alternatively been associated with jets, string ten- 
sion, and Heisenberg uncertainty quantum limits but not 
with collective final-state space-momentum correlations. 
The dAu system is unique in that the transverse momen- 
tum dependence of the HBT radii in three dimensions 
has only been studied at STAR. This system-wide com- 
parison is clearly important in understanding the physics 
behind the space-momentum correlations present in each 
system. 



Figure (0} shows the M T dependence of the three- 
dimensional Bertsch-Pratt Gaussian pion radii as seen in 
pp collisions. A significant dependence is seen in 
suggestive of a boost-invariant expansion along the beam 
direction. In addition, a small, but non- negligible, depen- 
dence of the transverse radii is also seen. In Fig. © the 
same pp data from Fig. Q are plotted in conjunction 
with the AuAu and dAu Gaussian radii (note the change 
of scale between Fig. (@J and 0). Like in Figure (2j), the 
system scale naturally grows smaller as one goes from 
central AuAu reactions to dAu and pp collisions. 

Most interesting is Figure which shows the ratio of 
the AuAu and dAu radii to the pp radii. The ratios are 
strikingly flat, given the very different mechanisms pre- 
sumably driving these systems. The similar R\ ong depen- 
dence is perhaps not unexpected since all of these systems 
were analyzed at mid-rapidity, a nearly boost-invariant 
region at STAR for all systems considered. Very different 
mechanisms, from hydrodynamic expansion to inside-out 
string fragmentation, give similar behavior of R\ ong ver- 
sus Mt strictly from boost-invariant arguments. 

However, in the transverse direction the results are 
more perplexing. In pp collisions, the transverse dynam- 
ics are most likely driven by the independent fragmenta- 
tion of (at most) a few strings, some of which will form 
jets. In AuAu collisions, the transverse dynamics are 
governed by collective expansion such as flow. It seems 
strange that these two mechanisms give rise to very simi- 
lar Mt dependence of the transverse radii. Perhaps much 
can be learned from this observation and the effect is still 
under investigation. 

There is a long history of using HBT in elemen- 
tary particle physics to study QCD and the space-time- 



7 ^STAR Preliminary 






STAR Preliminary 


7 


6^" 




Rout (fin) 








Rside (fin) = 


5 ■ 


■ 






■ 




^5 




a 






■ 


4^ 


■ 


■ 


■ 


■ 


■ 










■ 






1 :; 




■ 


■ 


■ 


■ 


■ 


■ ■ : 


f . 


T 

m 


▼ 

* ■ 


▼ 




V 

• 


-2 

t *t 










L-. , 1 , , , 




^0 



7 r. 

3^ 
2^ 
1r 



Rlong (fin) 



▼ 

■ 

STAR Preliminary 



Au-Au (central) 
Au-Au (midcentral) 
Au-Au (peripheral) 
d-Au 
PP 



0.25 0.3 0.35 0.4 0.45 0.5 0.55 

MT (GeV/c) 



FIG. 5: The transverse mass dependence of the 3D Bertsch- 
Pratt radii in pp, dAu, and AuAu at STAR. Three centralities 
are displayed for AuAu. The pp data points are the same as 
those shown in Figure (|p). 
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FIG. 6: The ratio of the dAu and AuAu radii to the pp radii 
as displayed in Figure (0. The lines are provided to guide 
the eye. 



momentum structure of hadronization. The multi-system 
capabilities of RHIC provide a unique laboratory for link- 
ing what has been historically studied in elementary par- 



ticle collisions and what is known about HBT in heavy 
ion physics. By studying the HBT of the pp and dAu sys- 
tems in the context of AuAu collisions at STAR, we hope 
to gain a better understanding of the freeze-out of nu- 
clear matter under various extreme conditions: from hot 
and dilute pp collisions, where the space-time- momentum 
structure of hadronization itself is probed - to the cooler 
and denser, highly interacting, nuclear medium gener- 
ated in AuAu reactions, where final state collective ef- 
fects dominate. 
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